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Exact positive and negative energy solutions for the eigenvalue problem of the Schrodinger equation 
in one dimension with a 5' interaction are found and analyzed. An infinite series of transparency 
resonance levels in the strength of this interaction is shown to exist. This result is against the 
actual belief that the 5' potential acts as a totally refiecting wall. A finite number of bound states 
is obtained, contrary to the previous result on the existence of only one bound state. A new efi'ect 
of a negative stepwise drop in the electron density across the 5' junction is observed. The solutions 
are also applied to the propagation of the electromagnetic field in dielectric media. 



PACS numbers: 03.65.-w; ll.lO.Gh; 73.22.Dj; 41.20.-q 
Point or contact interactions are widely used in various 
applications to quantum physics producing exactly 
solvable models of complicated physical phenomena . 
Intuitively, these interactions are understood as sharply 
localized potentials, exhibiting a number of interesting 
and intriguing features. Applications of these models 
to solid-state physics are of particular interest nowa- 
days, mainly because of the rapid progress in fabricating 
nanoscale quantum devices. Other applications arise, for 
instance, in dielectric media where electromagnetic waves 
scatter at boundaries or thin layers In this Letter, we 
report on a resonance transparency and a bound states 
structure of the derivative of Dirac's delta function poten- 
tial. This so-called 5' interaction has attracted consider- 
able attention because of a series of controversial results 
published in both physics- and mathematics-oriented lit- 
erature fl-[lOt|, many of which are not yet fully under- 
stood. In one dimension, and in the limit that neglects 
the interactions between electrons, the 5' interaction is 
formally given by the stationary Schrodinger equation 



^"{x) + a'^5'{x)il){x) = pV(a;), 



(1) 



where the prime denotes the differentiation over the spa- 
tial coordinate x, V'(a^) is the wave function for a particle 
of mass m (the units in which fi^ /2m = 1 are used), a 
the dimensionless strength interaction parameter, and 
(positive, zero, or negative) energy. 

It has been proposed ealier jj] that the spectral prob- 
lem of the Schrodinger equation with a point interac- 
tion at X = can be analyzed in terms of the equation 
{CqiI}){x) — p^iplx), Co = —dP/dx"^, for a; ^ 0, with ap- 
propriate boundary conditions at x = 0. As a result, a 
number of studies were devoted to finding correct bound- 
ary conditions at x = 0, where the derivative 6'{x) was 
located. The most striking feature of these results is that 
in some cases the boundary conditions appear to be ir- 
relevant to explicit solutions and one of these cases was 
analyzed in Ref. [Q. 

Since 5'{x) is a generalized function, it can be realized 
as a limit of regularized potentials, which we denote by 



Vi(x), with I being a regularization parameter, that con- 
verges to 5'{x) as Z ^ in the sense of distributions [ pT| . 
Then, solving the regularized version of Eq. , 

{CiA){x)=p'Mx), Ci = -d^ldx^ +a''Vi{x), (2) 

we analyze the behavior of the wave function (x) in the 
limit / — !■ 0. In this way, the subtleties of the boundary 
conditions at x = for the operator Cq can be avoided. 
Instead, a rigorous mathematical meaning to the limit 
of the sequence {C{\ in the Hilbert space i2(— oo,oo) 
must be given. To this end, the notion of the graph 
limit introduced ealier by Glimm and Jaffe to prove 
the self-adjointness of the Hamiltonian for the relativistic 
Yukawa interaction in one spatial dimension, can be used. 
Thus, let {V'i} be a sequence of vectors such that there 
exist "0 — linii-^o "0; ^-nd 9 = limi^o('CiV'i)j forming a 
pair of limit vectors {0, 0}. Since the domain of each Ci 
is dense in L2(— oo,cx)), the limiting operator C called 
the graph limit of {Ci\ is defined by putting = £0. 

In our special case, it is convenient to construct a se- 
quence of regularizing functions {V;(x)} using step ap- 
proximations to the 5' {x) function. First, based on phys- 
ical intuition, we are able to decide how close is the 6' 
potential to the actual (regular) potential. Second, our 
procedure demonstrates how to realize explicitly the wave 
function discontinuity at x = 0, the problem of recent 
interest Third, step- like potentials can easily be 

manufactured using, e.g., thin layers of different types of 
semiconductors. More specifically, the 5' interaction may 
be used as a limiting case of a device in which a small 
region of large repulsive potential is followed by a small 
region of large attractive potential, which may be called 
a barrier-well junction. The similar experimental situa- 
tion can be arranged in dielectric media Therefore, 
the results reported here are applicable to larger classes 
of wave equations than the Schrodinger equation. 

It is instructive to consider two types of regularizing se- 
quences {V;(x)}, which are expected to result in different 
physics in the limit I 0. The key point is to impose dif- 
ferent constraints on the behavior of {Vi(x)} in the vicin- 
ity of X = 0. For the first type, each function [denoted by 
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is supposed to be identically zero in the vicinity 
of a; = 0, putting, e.g., V,°(a;) = [S{x + I) - 6{x - l)]/{2l) 
0, or appropriate step functions, the support of which 
does not contain the point x = Q. Then for each / > 0, 
we have an oscillatory solution in the vicinity of cc = 0, 
which collapses to a continuous function with a node at 
a; = when / — > 0, so that an incident current is totally 
reflected at this point 0. Note that for each / > 0, the 
energy flow is partially transmitted. For the second type 
of regularizing functions Vi'^{x), an opposite constraint 
at a; = 0, namely the existence of a discontinuity at this 
point that goes to infinity as I ^ 0, is imposed. As an ex- 
ample, the step function VJ''(x) shown in Fig. |l| that splits 
the X axis into four regions: Vi'^{x) = (|a;| > Z > 0), l'^ 
{—I < X <0), and — {0 < x <l), can be chosen. With 
this potential, at the x = point, two different types of a 
solution are connected, namely tunneling and oscillatory 
ones. Therefore for each finite /, the wave function ipi{x) 
has no node at x = 0, and one could expect that at least 
for some system parameter values, the limit of continu- 
ous functions 'ipi {x) will have a finite discontinuity at this 
point. This would immediately implies a nonzero current 
flow across the 5' junction as shown below. 

The current (energy flow) transmission coefficient T is 
calculated from the scattering solution ipii^) of Eq. (||) 
with > in a standard way. For the potential depicted 
in Fig. |l|, it appears to depend on the two dimensionless 
parameters rj = Ip and a: 

T = 4 |2 + [cos(a) cosh(/3) - (a//3) sin(a) sinh(/3)]^ 

-t- [cos(a) cosh(/3) + {(3/a) sin(Q:) sinh(/3)]^ 

-I- 77^ [a~^ sin(a) cosh(/3) -I- /3"^ cos(a) sinh(/3)] ^ 

-|- r/^^ [asin(a) cosh(/3) — /3cos(a) sinh(/3)]^| , (3) 

where a = cr"^ + rj'^ and f3 = \J — rf' . Note that the 
transmission across this junction does not depend on the 
sign in front of in Eq. ([|), i.e., it is the same both 
from the left and from the right. As shown in Fig. ||, 
the transmission over the barrier with energies 77 close to 
the barrier suffers big variations while the parameter u 
changes. 

Let us consider the limit 77 — > that corresponds to 
the limit I 0. In this case, a,/3 cr, while T 0, 
except for those values of a that satisfy the equation 

tan a = tanh a. (4) 

This equation, except for the trivial solution ctq = 
admits a discrete countable set of "resonance" levels 
cr„, n = 1, 2, . . . , at which the 8' junction becomes par- 
tially transparent, in contrast to the general belief that 
the 6' interaction acts as a reflecting wall j?). Indeed, in 
the limit — > 0, from Eq. (||) one obtains the limiting 
values of the transmission coefficient T: 



Tn = lim T{r], cr„) = 1 — tanh^ (t„. n = 1, 2, . . . . (5) 

77^0 

One can see from the last equation that there is no de- 
pendence of T on the dimensionless energy parameter 
77. Its behavior is illustrated in Fig. || as a function of 
the amplitude parameter a for different values of the di- 
mensionless parameter rj. It is clearly shown that the 
continuous curve T — T{ct) as a function of a tends to 
isolated, descreasing (in a) peaks as 7; ^ (see Fig. ||). 

One can easily find the left and the right values of 
the limiting wave function ^{x) as x — > ±0. Let us de- 
note these values of the wave function ili{x) from the 
barrier and the well sides as ij^b = lim/^oV'/(^0 
V'u, = lim;^o'0;(O: respectively. 

Using the scattering solution of Eq. (||) with > 0, 
one finds the ratio 

\^b/ipw\^ = (1 - tanh^ an) /{I + tanh^ cr„) < 1, (6) 

which shows that the limiting wave function rjj{x) has 
a finite nonzero jump at x = if cr = (t„. The last 
inequality means that the probability to find an electron 
nearby the junction is higher from the side of the well 
than that of the barrier. In other words, at the resonant 
levels {(Jn}^=i, a negative drop in the electron density 
occurs, while passing across the junction from the barrier 
to the well side. 

In order to find a discrete spectrum of bound states, 
we solve Eq. (^) with negative energy p^ < 0, putting 

= ly^—p"^. As a result, one finds the equation on the 
discrete spectrum C„, 71 = 0, 1, ... : 

(C^/a/3) tan(Q!) tanh(/3) + (C/2) (a~HanQ; + (3'^ tanh/3) 
+ 1 - (atana - /3tanh/3)/2C = 0, (7) 

where a — \/ o"^ ~ and /3 — \J + Cp' ■ The solutions 
of Eq. (|^) with respect to ( depend on a and they form a 
finite discrete spectrum {Cn}n=o depicted in Fig. ^. The 
number of bound states N ~ N{a) is determined by the 
inequalities aN < a < CAf+i. In the vicinity of small C, 
Eq. (0) is asymptotically reduced to 

a tan a — /3 tanh /? = 2^, (8) 

from which one finds that in the limit ^ ^ (Z — > 0), the 
77,th eigenvalue C,n{o'), = 0, 1, . . . , begins at ct = (t„, i.e., 
C,n{(^n) = as shown in Fig. |^. Note that each eigenfunc- 
tion ipnix) that corresponds to the nth. eigenvalue with 
a = On has at a: = the same discontinuity jump as the 
scattering wave function ^{x) for the same value cr„. 

In quantum mechanics applications, the dimensionless 
quantity tr that scales the interaction strength is a sys- 
tem parameter and therefore cannot be controlled during 
electron transport across the 5' junction. However, in 
integrated optics , this parameter can appear as a fre- 
quency of the electromagnetic field. Moreover, the results 
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of the spectral analysis obtained above for the 5' interac- 
tion can be applied there directly. Thus, let us consider 
the Maxwell equations in a medium with one-dimensional 
(say, along the x axis as shown in Fig. ^ inhomogeneity 
of dielectric permeability. In the Cartesian coordinate 
system, for TE waves {E^ = Hy = = 0) with fre- 
quency CO, the Ey component satisfies the equation 



d^E„ d^E, 



dz 



}■ + k^e{x)Ey = 0, 



(9) 



where k — uj^eofi is the wave number with Eq being 
the dielectric permeability of vacuum and /i the mag- 
netic permeability, and e{x) is a dimensionless profile 
of dielectric permeability along the x axis. Assuming 
next a wave propagation along the z axis, i.e., putting 
Ey{x,z) — E(x) expiiqz) where g is a wave vector, one 
obtains the eigenvalue problem: 



E"{x) + k^e{x)E{x) = q^Eix), 



(10) 



with two spectral parameters: wave numbers k and q. 

We consider the e{x) profile of the form shown in Fig. ||, 
where e{x) = if |a;| > a, and e{x) = eb+£m if < a; < a 
and £{x) = £b — £m if —a < x < 0. Using the scaling 
transformation x ax/l and E{x) — > E{ax/l) = 'ijji{x), 
Eq. ( p^ ) can be reduced to Eqs. (||), in which 



k^a Ern and ^ {a/iy{k^eb - q ). 



(11) 



Then in Eqs. (|^) and (|7|), we have to put j] — a\J k^Sb — q^ 
and C = a\J(p — W^b , respectively. Therefore, the 
cut-off frequencies lu (or k) that cut the discrete spec- 
trum of waveguide modes from the continuous spectrum 
(77 = C = 0) lie on the line q — k^feb-, as shown by line 
1 in Fig. ^ More precisely, they are given by the points 
(fcni'Zn) on the (fc,g) plane, where fc„ = On/ [a^/e^) and 
In = {'^n/a)y'eb/£m , with (T„'s, n = 1, 2, be- 
ing positive roots of Eq. (^). Moreover, according to 
Eq. (^, one can conclude that the cut-off and the trans- 
parency frequencies coincide. In other words, the junc- 
tion of the form shown in Fig. ^ is partially transparent 
in the transverse direction at the cut-off frequencies of 
the waveguide regime that occurs in the longitudinal di- 
rection. These results are illustrated by Fig. Using 
Eqs. (0), one finds that the waveguide regime is given 
by the set of curves (?„(fc) = y/k'^Sb + Cn (ka^/e^) / 



which are depicted above line 1 in Fig 
splits (see line 2 in Fig. 



\ The equation 
P the whole sec- 



q = ky/Sb - e _ 
tor of continuum spectrum into two sectors I and II that 
correspond to 77 < cr and rj > a, respectively. In these 
sectors, the transparency properties given by the trans- 
mission coefficient r(ry, a) are demonstrated in Fig. |^. 

Thus, using the sequence of physically motivated reg- 
ularizing functions shown in Fig. ^, an infinite, count- 
able set of resonance values for the interaction strength 
parameter a, at which a transparent regime of current 



flow occurs, has been found explicitly for the 6' inter- 
action described by the Schrodinger equation (|l|). All 
these resonance levels, {cr„}J^i, appear to be positive 
roots of the simple equation (Q). This result contradicts 
Patil's result , according to which the transparency is 
identically zero. The reason of this controversy emerges 
from the fact that the approximating sequence {VJ°(a;)} 
constructed from the double of Dirac's S functions 0, 
necessary results in a continuous wave function with a 
node at a; = as ^ — *■ 0, blocking a current across this 
point. A controversy also arises concerning the bound 
states structure. Instead of one bound state as claimed 
previously 0, a finite number of such states is found, 
increasing with a. The profiles of the bound states are 
shown to have also finite jumps at x = 0. 

Finally, note that the graph limit C cannot be rep- 
resented as the sum of Cq and a well-defined object 
S'{x)il){x). Since ip{x) is a discontinuous function at 
X — 0, this product is meaningless, but in the regulariza- 
tion process, the infinite renormalization of Vi'''{x)ipi{x) 
occurs, so that divergent terms in the limiting operator 
C are totally cancelled. 
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I. FIGURE CAPTIONS 



FIG. 1. The piecewise constant approximation of the S'{x) 
function, where I is the width (regularization) parameter. 

FIG. 2. Dependence of the transmission coefficient T on rj 
and a: (a) Three-dimensional surface T{ri, a), where the white 
curve shows the values of T on the line rj = a. (b) Section 
of the surface T{rj, a) a.t rj — 15. The inset demonstrates big 
variations of the coefficient T nearby = a. 

FIG. 3. Logarithm of the transmission coefficient T 
against a at rj — 0.0005 (solid curve), 0.05 (dashed curve), 
and 0.5 (dashed-dotted curve). 

FIG. 4. Five bound states n — 0,1, 2, 3, 4, as functions 
of a depicted up to cr = 15. 

FIG. 5. Profile of the dimensionless step dielectric perme- 
ability e{x) along the x axis; 2a is the width of the junction 
layer, et a background dielectric permeability, and Sm. the 
amplitude of permeability variation across the layer. 

FIG. 6. Discrete and continuous plane spectra of electro- 
magnetic wave transmission across the layer structure shown 
in Fig. ^ The cut-off frequencies {k„,q'^) lie on the line 
q = k^/Fb labeled by 1. Above this line, the first four waveg- 
uide modes are shown by curves qnik), n = 0, 1,2,3. Below, 
in the continuous spectrum sector, the line q = k^Jsb — £m 
labeled by 2 represents the scattering regime nearby the min- 
imum of dielectric permeability, e — Sh Em • 
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This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/math-ph/0204045v 1 
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